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Quantum Monte Carlo is used to calculate various pairing correlations of the 2D Hubbard model 
possessing band features experimentally observed in the cuprates. In the hole-doped case, where 
the Fermi level lies close to the van Hove singularities around (0, 7r), the d x 2_ y 2 pairing correlation 
is selectively enhanced, while in the electron-doped case, where the singularities are far below the 
Fermi level and the Fermi surface runs through (±7r/2, ±7r/2), both d x 2_ y 2 and d xy correlations are 
enhanced with the latter having a y/2 x y/2 structure. The two pairing symmetries can mix to result 
in a nodeless gap. 

PACS numbers: 74.20.Mn, 71.10.Fd 
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Since the seminal proposal by Anderson, jj| great the- 
oretical effort has been made to investigate the possibil- 
ity of describing various aspects of the high Tc cuprates 
within the two dimensional (2D) Hubbard model. ||] 
Those include the antiferromagnetic insulating phase in 
the undoped systems, the normal state above Tc, and 
the superconducting state. Among all, it has been an 
open question whether the Hubbard model can actually 
account for the superconductivity, especially its pairing 
symmetry, in both hole and electron doped cuprates. 

Many analytical calculations have supported the pos- 
sibility of d x 2_ y 2 pairing in the nearly half- filled 2D Hub- 
bard model. ||-|(| While some previous numerical studies 
of pairing correlations in finite systems have given nega- 
tive results for superconductivity, we have recently 
shown that an enhanced d x 2_ y 2 pairing correlation is in- 
deed detected numerically if we ensure that the highest 
occupied one-electron levels (HOL) and lowest unoccu- 
pied levels (LUL) at U = in finite systems are suffi- 
ciently close. [0 This precaution, as motivated from the 
numerical studies on Hubbard ladders, Jll],[l^] has been 
necessitated because the energy scale of the superconduc- 
tivity in the Hubbard model, if any, should be of the order 
of O.Oli, HU while the discreteness of the energy levels in 
finite systems tractable in numerical calculation is much 
larger (~ O.lt) unless parameter values are tuned. 

In the hole-doped cuprates such as YBCO and BSCCO 
there is now a body of accumulating evidence that the 
pairing symmetry is d x 2_ y 2 (at least around the opti- 
mal doping), Jd| which is consistent with the previous 
Hubbard-model studies. On the other hand, experimen- 
tal results for the electron-doped NCCO seems to indi- 
cate an s-wave, or more precisely, a symmetry with a 
nodeless superconducting gap. ]T^-p^| Experiments have 
also revealed further differences between hole-doped and 
electron-doped systems. Specifically, the angle-resolved 
photoemission spectroscopy (ARPES) has shown that 
the 'extended' van Hove singularity (VHS) around k = 
(0, 7r) and (it, 0) lies very close to the Fermi level (Ep) 



in YBCO and BSCCO, while the VHS lies far below the 
Fermi level in NCCO. @ 

The purpose of the present study is to explore whether 
the difference of the pairing symmetry between electron 
and hole doped systems can be explained within the 2D 
Hubbard model possessing the band features observed 
experimentally. The essential band features (namely, the 
shape of the Fermi surface and the relative position of the 
VHS to E F ) of YBCO, BSCCO, and NCCO can be re- 
produced by introducing a next-nearest neighbor (NNN) 
transfer about half the nearest neighbor (NN) one. Note 
that in our previous study mentioned above [ ^0[ , such 
an electron-hole asymmetry was not taken into account 
since we considered the Hubbard model with only NN 
transfers. 

Quantum Monte Carlo (QMC) method is used to cal- 
culate correlation functions of d xy , NN and NNN ex- 
tended s as well as d x 2_ y 2 pairings. To look into such 
diverse symmetries has been motivated from the follow- 
ing physical consideration. Namely, the pair-tunneling 
processes between (ki |, — ki |) and (k 2 j, — k 2 |) favors 
the pairing order parameter A that satisfies A(ki) = 
— A(k 2 ), which is a picture known to be at work in the 
two-leg and three-leg f|l]j2|-§|] Hubbard ladders. 

In this picture, the pair-tunneling between the fc-points 
around (0, w) and (tt, 0) favors d x 2_ y 2 pairing. Such pro- 
cesses should indeed be pronounced in the hole-doped 
cuprates because Ep lies close to (0, it) and (tt, 0), and 
the density of states around these points is large. Thus, 
in this case, d x 2_ y 2 pairing should be dominant with pos- 
sibly other symmetries mixing slightly. By contrast, in 
the electron doped case, other pair-tunneling processes 
may set in on a nearly equal footing in determining the 
pairing symmetry, because VHS lies far below Ep. Then, 
not only d x 2_ y 2 but also d xy pairing or extended s pair- 
ing (with gap functions that have nodes on the Fermi 
surface, but do not change sign by a 90 degree rotation) 
will become eligible, so that some of these symmetries 
may mix with comparable weights. 
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In fact, we find here that the d x 2_ y 2 correlation is dom- 
inant in the hole-doped case, while in the electron-doped 
case both d x 2_ y 2 and d xy correlations are enhanced, with 
the latter having a v2 x \[2 structure. These two pair- 
ing symmetries can in fact mix ending up with a nodeless 
gap without breaking the time reversal symmetry, unlike 
in d x 2_ y 2 + id xy pairing p6| , p7t where the symmetry is 
broken. Correlation of the extended s-wave pairings is 
found to be suppressed in all the cases investigated. 

We consider the 2D Hubbard model on a square lattice 
with NN (t), NNN (t'), and third NN (i") hoppings, 

^ = — y ] \j'^( c x.y,cr c x + l,y,cr) + ty (c x ^ y .o-Cx,y+l,o-) 
x,y,cr 

+ t'/{ c x,y,a c x+l,y+l,(r) + t'\{ c x,y,a c x-l,y+l,a) 
+ t "( c l,y,a c x+2,y,a- + cl^ y ^C x . y+2 .a) + h.C.] 

x,y 

Here, (x, y) is the coordinate of the sites, and the lattice 
constant is taken as unity. Periodic boundary condition 
is assumed, and we set t r = 1 hereafter. 
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As mentioned above |10|, it is necessary to put Ep 
at U = between the HOL's and LUL's separated by 
an energy of Ae° less than 0(0.01) in order to detect 
a symptom of superconductivity having an energy scale 
of O{100K). On the other hand, QMC is unstable for 
exactly As = 0, namely for open shell configurations. 
Thus, we accomplish Ae a ~ 0(0.01) by making t x and 
t y , and/or t', and i( slightly different, where t x ^ t y lifts 
the degeneracy between (±ki,±k2) and (±/c 2 ,±&i) for 
I fci| 7^ 1^1, while t' i ^ t'^ lifts the degeneracy between 
(fcijfci) and (±fci,=Ffci). 

We have employed the ground-state, canonical- 
ensemble QMC, where we have implemented the stabi- 
lization algorithm adopted by several authors. J2f| We 
adopt the free Fermi sea as the trial state, and take the 
projection imaginary time r up to ~ 40 to ensure the 
convergence. Small Ae° makes the negative sign problem 
serious, but by taking a relatively small value of U(= 1), 
we can check the convergence with respect to r without 
running into serious sign problem. 

We have calculated the pairing correlation functions, 

P ( r ) = (0\x + Ax,y + Ay)0(x,y)} with 

\Ax\ + \Ay\=r 

Onn(£, V) = 2_j a ( c x,y,crC x +S,y,-a- ± Cx,y,<rCx,y+S,-a) 
5,a 

OnNn(^,2/) = ^ &(c x , y ,crC x+ s, y +S,-cr ± Cx,y,<TCx—S,y+S,-<r)) 

where 8 = ±1. The plus (minus) sign in Onn corresponds 
to NN s (d x 2_ y 2) symmetries, while the plus (minus) in 
Onnn to NNN s (d xy ) symmetries. 
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FIG. 1. QMC result for d 
correlations for a hole-doped system 
118 electrons with n = 0.82). t y = 



(b) pairing 
12 x 12 system with 
0.999, t'. 



I 



-0.429, 



t[ = -0.43, t" = 0.07, and U = 1 (O)- The dashed lines 
represent the U — result. The inset shows the HOL's and 
LUL's within 0.01 to E F . 



We have looked into various values of n, t', and t" in- 
cluding other than the ones described below, and found 
NN and NNN s-wave pairing correlations to be strongly 
suppressed with U at large distances. At first this may 
seem odd because these pairings do not have any on-site 
amplitude. This might be because the extended s pair- 
ings always couple, at least at the mean-field level, with 
the on-site s pairing, [ p9[ which is directly suppressed 
with U > 0. Thus, we show only d x 2_ y 2 and d xy pairing 
correlations in the following. 

We first look at the hole-doped case. We consider 
a 12 x 12 lattice with 118 electrons (band filling n = 



-0.429, t{ 



0.82) with t y = 0.999, t', 
t" = 0.07. For this choice, the HOL's 



12) 



= -0.43, and 
at (±7r/6, 7t), 
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(±7r/6, ±5tt/6), (7r/3,-7r/3), (-7r/3,7r/3) and LUL's at 
(tt, ±7r/6), (±57r/6,±7r/6), (tt/3, tt/3), (-tt/3, -tt/3) lie 
within 0.01 in energy at U = 0. The Fermi surface, repre- 
sented by these HOL's and LUL's is displayed as an inset 
in Fig.|l|(a). There, reflecting the high density of states 
around VHS, many /c-points around (0,7r) and (tt,0) ap- 
pear, while the points around \k x \ — \k y \, although fewer, 
also exist. 

In Fig.[l](a), we show the d x 2_ y 2 correlation as a func- 
tion of real space distance r = |Ax| + 1 Ay\. It can be seen 
that the correlation is enhanced for U = 1 over that for 
U = 0, especially at large distances. By contrast, the d xy 
correlation shown in (b) is not enhanced within the error 
bars. The dominant d x 2_ y 2 pairing is consistent with the 
expectation from the pair-tunneling picture given above. 
On the other hand, we cannot rule out the possibility of 
a small d xy mixing, since if more fc-points exist in the 
vicinity of Ep, not only the d x 2_ y 2 correlation would 
be more enhanced, but also the d xy might be enhanced, 
which would imply their mixture. Further, d xy may mix 
in a time-reversal broken form, d x 2_j ( 2+id ;r j ) , |2^,|2^] espe- 
cially in magnetic fields, (3^] which is of interest from the 
viewpoint of the recent experimental observations sug- 
gesting such a possibility at low temperatures. 32 33] 

Let us now turn to the case of electron doping. This 
time, we take 190 electrons /12 x 12 (n = 1.32) with 



0.999, t' 



-0.499, t[ 



-0.5, and t" 



0. (In the actual calculation we have employed the 
electron-hole transformation to consider a 98 electron 
system with t' > 0). Here, HOL's reside at (±7r/3, tt), 
(7r/2,-7r/2), (-71-/2,71-/2), while LUL's at (7r,±7r/3), 
(7r/2,7r/2), (-7r/2,-7r/2) for U = (inset of Fig.|(a)). 
Note that (±7r/2, ±7r/2) lies right on the Fermi surface, 
a feature seen in the ARPES data of NCCO. @ 

The QMC result in Fig.^(a) shows that, although the 
Fermi surface is now shifted away from (tt, 0), (0, tt) down 
to (71-, ±7r/3), (±7t/3,7t), we still have an enhancement 
of the d x 2_ y 2 correlation, although the enhancement is 
smaller than that in the holed-doped case. 

Now, more striking is the behavior of the d xy corre- 



xy 



lation shown in Fig. 2(b). At large distances, the d 
correlation is enhanced at even distances (Ax + Ay = 
even), while suppressed at odd distances, which means 
that it has a \f2 x \/2 superstructure. A Fourier trans- 
form of the correlation function indeed shows that its 
(71-, tt) component is enhanced with U. 

The result suggests a coexistence of the d x 2_ y 2 and 
the \J~2 x y^2 d xy pairings, whose order parameters are 
CkT c -ki(cosfc x - cosk y ) and c k TC_( k +Q)x(sin k x sin k y ), 
respectively, with Q = (tt, tt). If they both have 
long-range orders, we should take (ck-f,Ck+Qt) an d 
(c-ki, C-(k+Q)i) as basis to diagonalize the 2x2 order 
parameter matrix to have 

A±(k) = ±y^4(cos k x — cos k y ) 2 + B(smk x s'mk y ) 2 , 



where A, B > 0. This form, which is nodeless, is similar 
to the energy spectrum of the chiral spin state proposed 
by Wen, Wilczek, and Zee. Q The order parameter of 
the chiral spin state is defined for (c^c), the hopping am- 
plitude, while we are here talking about (cc), the pairing 
amplitude. The corresponding superconducting gap co- 
incides with that of the d x 2_ y 2 +id xy pairing, p6| , p7t but 
we must stress that the present order parameter is real 
and hence does not break the time reversal symmetry as 
in d x 2_ y 2 + id xy . Thus we end up with a fully-gapped, 
time-reversal-symmetric mixture of d x 2_ y 2 and d^j, pair- 
ings. 

As seen in Figj| (b), the V2 x V2 structure of the d xy 
correlation is not observed in the hole-doped case. In 
fact, we have considered a wide variety of cases, some of 
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FIG. 2. A plot similar to Fig.[j] for an electron-doped sys- 
tem (190 electrons /12 x 12 with n = 1.32) for t y = 0.999, 
-0.5, t" = 0, and U = 1. 



t'/ = -0.499, t[ 
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which will be published elsewhere, and found that the 
structure in the d xy pairing emerges only when 
(±7r/2, ±7r/2) lies on the Fermi surface. Then, the dif- 
ference in the pairing symmetry between the hole-doped 
and electron-doped cases may be not only due to the 
relative position of the VHS against Ep, but may also 
come from the fact that (±7r/2, ±7r/2) lies very close to 
the Fermi surface in NCCO. 

The relation of (±7r/2, ±7r/2) to the pairing having a 
superstructure has also been suggested for the t-J model 
by Ogata quite recently Q. Using a variational ap- 
proach to the t-J model, he showed that the energy of 
& x 2_ y 2 pairing state is lowered with a full gap when mixed 
with KK cxtcnded-s pairing having finite momentum of 
(it, 0) or (0,7r), if the system is lightly doped, so that 
(±7r/2,±7r/2) is close to the Fermi surface. §7| The dif- 
ference with the present d xy superstructure is that the 
pair here has a finite total crystal momentum of (tt, tt) 
resulting in a \J~2 x \/2 structure, while Ogata's extended 
s-pair has a momentum (tt, 0) or (0, tt) with a 2 x 2 struc- 
ture. As mentioned above, we have so far found extended 
s-wave correlations to be suppressed at large distances, 
but we believe further calculation for various values of 
parameters is necessary to reveal the relation between 
the present result for the Hubbard model and Ogata's 
result for the t-J model. 

In summary, we have shown that the 2D Hubbard 
model possessing band features experimentally observed 
in the cuprates can account for both the d x 2_ y 2 pair- 
ing for hole doping and a nodeless pairing for electron 
doping. The fact that the present result is obtained for 
rather small values of t) suggests that large inter- 
actions (U ^> t) may not be essential to the occurrence 
of superconductivity, although the strength of the inter- 
action will certainly dominate the absolute magnitude of 
the gap or Tp. 
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